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Overview 


- Volume representations 
■ Voxelization survey 

- Lessons from Windborne 


mi 


Who has worked on a voxel engine before? Who wants to? 

My goal is to give the talk I wish I would have had before I 
started on our procedural engine. 


Three parts to this talk. A lot of content, so I'll move fast, but 
my goal is to give a background on what to look up (that was 
the biggest problem I ran into when first approaching this 
subject). You have to KNOW that you should look up contour 
indicator functions for instance. 
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Applications 



Image credit: 

http://www.pcgamer.com/everquest-next-landmark-hands-on/ 



Image credit: 

http://www.minecraftforum.net/news/60286-original-minecraft- 

surpasses-22-million-sales 



Image credit: 

http://pixologlc.com/zbrush/gallerY/ 


Image credit: 

http://www.autodesk.com/products/3ds-max/ 


Image credit: 

Defense Grid level editor (Hidden Path) 




Procedural geometry is useful for far more than just dynamic 
worlds as people first think- really all 3d art creation requires 
it. 

While we normally don't deal with this as engineers, knowing 
the methods and having a rich toolset can help us make better 
games. 
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Part 1 

Volume Representations 


LI ■ ■ 


4 


GDC 




GAME DEVELOPERS CONFERENCE March 14-18, 2016 Expo: March 16-18, 2016 0GDC16 



Quaternions, a love story 




Smiley Cube, 
random Quaternion 
transform 


Smiley Cube, 
random Matrix 
transform 


ill 


Everyone knows Quaternions are good for working with 
rotations, and Matrices have problems. 


The issue I want to call attention to is that matrices are 
overspecified. There are lots of things that matrices can 
represent that aren't normalized orthogonal basis sets, 
whereas that's all quats can represent. 

If I fill a quaternion randomly with data, it's good. If I fill a 
matrix randomly, it's not good. 
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Surface Representations 



Manifold Non-Manifold Non-Manifold 

Image credit: 

http://www.autodesk.com/techpubs/allasstudlo/20l0/index.htm 
l?uri=WS73099ccl42f487557230b50811d7dl92c64- 
79bl.htm, topicNumber=d0e884l9 


Orientation 

Does c lie on. to the left of, 
or to the right of aB'l 



a x a \ 

b x b y 

c x C v 


\o x -c x a -c l 

\ b '~ c , *v" c J 


Incircle 

Does d lie on, inside, or 
or outside of abcl 

a,-d x Oy-dy (a,-«4) 2+ (Ot- 4^1 

■ h t -d x by-dy 

c t -d, Cy-dy (c,-4,) 2 +(<;»-4r) 2 | 


Image credit: 

https://www.cs.cmu.edu/~quake/robust.html 


a, a, «,'+ I 

b, ft, b?+by 2 I 

Cy Cy Cyb Cy I 

d x d dc+dy 2 1 



Surface representations have the same problem matrices do. 
They are subject to euler's identity and manifold issues, 
exacerbated by floating point error. 

It's very hard to make a left right plane test that's 100% 
robust to floating point error. You WILL hit problems if you roll 
your own surface rep CSG. 


Validation + Repair can be very complicated 
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Volume Representations 




■ Any continuous f(x,y,z) 

■ Implicitly define volume as f(x,y,z) > 0 

■ Surface is level set where f(x,y,z) = 0 



F(x,y) = Sin(x A 2)+cos(y A 2)-l Level set (c=0) 

Image credit: wolfram alpha 


i ■ n 

Function can be ANYTHING, no robustness problems- it's 
strictly a sampling problem. No verify and fix step. 
Programming without worrying about edge cases, just provide 
a function. 

Each point in space has a very concrete 'in' or 'out' and it's 
based on a continuous function, so it's manifold. 
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Volume Representations 




- Just need a continuous function 


■ Any algebraic function 

■ Signed distance field 


■ CSG trees (http://www.merl.com/publications/docs/TR2006-054.Ddf) 


■ Trilinear interpolation over a grid 

■ Density function 

■ Trilinear interpolation over a grid 



Image credit: 

https://en.wikipedia.org/wiki/TrillnearJnterpolation 


Or more complicated interpolation works 
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Density vs Signed Distance 


Density is much easier to work with (changes are local), but 
you lose out on some useful operations on the distance field 
(ability to expand/shrink volume and higher quality gradient 
calculations) 

You can treat density as a distance field with a clamped 
distance of approximately one sampling grid cell. 

Note: the level set / surface for both these methods is almost 
identical (mid grey in the bitmaps) 


Density 


Distance Field 



Image credit: http://forum.unity3d.com/threads/unity-signed- 
dlstance-field-textures.61243/ 


Image credit: http://forum.unity3d.com/threads/unity-signed- 
distance-fleld-textures.61243/ 


Link 
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Rendering 


Volume Rendering 



Image credit: 

https://en.wiklpedta.org/wikl/Ray_tracmg_(graphlcs) 



Convert to Surface 



Image credit: 

http://www.cs.carleton.edu/cs_comps/0405/shape/mar 

ching_cubes.html 


L II 


Ray tracing - evaluate function along ray, looking for 0 
crossings 

Convert to surface is what I'm going to talk more about 
because it fits better in the full volume of video game 
rendering knowledge we have now 
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Voxelization 

1. Sample f on a grid 

2 . Approximate the 
surface in each cell 

3. Make sure surfaces 
align at cell boundaries 

4. Profit! 



We used density fields because they are more dynamically 
updatable 
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Part 2 

Voxelization Survey 
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Voxelization Ideals 

■ Easy to Implement 

■ Local Independence 

■ Smooth 

■ Adaptive / suited to LOD 

■ Minimize triangle slivers 

■ Preserve sharp features 

■ Preserve thin features 





Image credit: 

http://aptww.org/IntlCatalog.nsf/vPlaylistItemByID/703 
45- 3/$FlLE/tree large. Jpg 


L II 


Low intercell dependencies if we want to be dynamic 
(remember, density fields over distance fields) 

Reduced triangle count - ideally large flat areas aren't an 
array of polygons 

Preserve sharp features (trilinear filter is a blur). We can't 
represent anything over the nyquist frequency of two grid 
cells. 



Simple Cubes 


■ Sample f(x,y,z) in the middle of each grid cell 

■ Draw a face between grid cells with different signs 
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The nice thing is we only have polies on the surface, so it's 
pretty efficient 
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Marching Cubes 

http://www.eecs.berkeley.edu/~jrs/meshpapers/LorensenCline.pdf 


■ Sample f(x,y,z) in the corner of each grid cell 

■ Use signs of corners for topology of triangles 


■ Locate verts at interpolated 0 point along edges 







Y 





yX— 

/ 










Image credit: 

http://www.eecs.berkeley.edu/~jrs/meshpapers 

/LorensenCline.pdf 


If you only look at the sign of the corners of each grid cell, 
there are only 15 topologies to choose from. You can slide the 
points along the edges based on relative weights if you have 
more than -1/1 
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Marching Cubes 


Easy to Implement: + 
Local Independence: 4- 
Smooth: + 

Adaptive: - 
Minimize slivers: - 
Sharp features: - 
Thin features : - 


Image credit: https://cl.iggcdn.com/lndlegogo-media- 
prod- 

cld/image/upload/c_fill,f_auto,h_4l3,w_620/vl406808 
220/ vv4nvjafctgzemwbevua.jpg 
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Transvoxel Algorithm 

http://transvoxel.org/ 



♦ 

m *** Y h S * * 


■ Sample f(x,y,z) in the corner of each 
grid cell 

■ Allow subdividing grid cell sides once 
(to seam adjacent LOD levels) 

• Use signs of sample points for 
topology of triangles 

■ Locate verts at interpolated 0 point 
along edges 


JO. tsi $ 1*3 M Mf iff jmt 

V* tmf yh « ty * *> f 

Jsr. Ja’ Va ^ 

-iff 'M- iff iff yS : 

■iff g g g] 

•jW £ ijf £ '• 

y* > -iff' £» W Jfef fr * 

**;: ?L' g 4,' 

Image credit: http://transvoxel.org/ 


Transvoxel is a method to allow marching cubes to span 
different LOD levels. 71 total topologies to handle a 
tessellation of any combination of sides 


GDC 


<> 


GAME DEVELOPERS CONFERENCE March 14-18, 2016 Expo: March 16-18, 2016 0GDC16 


Transvoxel Algorithm 




Easy to Implement: + 
Local Independence: 4- 
Smooth: + 

Adaptive: + 

Minimize slivers: - 
Sharp features: - 
Thin features : - 


Image credit: http://the31stgame.com/ 
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Dual Contouring 

http://www.frankpetterson.com/publications/dualcontour/dualcontour.pdf 


■ Sample f(x,y,z) in the corner of each grid cell 

■ Sample f'(x,y,z) at the location of each edge intersection 

■ Find an 'ideal' point within each cell, on the contour 

• Connect dual points for neighboring grid cells (supports multiple 
LOD resolutions) 



Dual contouring is a method to both spawn LOD levels and 
preserve sharp features. 

It's a dual method, so instead of create vertices on the edges 
of the grid, we create vertices inside the cells and connect 
them 
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Dual Contouring 



Easy to Implement: - 
Local Independence: - 
Smooth: 4- 
Adaptive: + 

Minimize slivers: + 
Sharp features: + 

Thin features : ™ 


Image credit: http://www.voxelfarm.com/ 
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Dual Marching Cubes 

https://www.cs.rice.edu/~jwarren/papers/dmc.pdf 


■ Sample f(x,y,z) over a fine grid 

■ Find the point that minimizes error (QEF) 

■ Subdivide octree at that point if error > £ 

■ Repeat 1-3 until error < £ everywhere 

■ Construct topological dual of this octree 

■ Tessellate (as Dual Contouring) 




Marching Cubes 


Dual Marching 
Cubes 


This extension realizes that marching cubes can be ran on the 
dual of an octree itself. Fit the grid to the data. 

Minimizing slivers is a proposed extension. 



Dual Marching Cubes 


Easy to Implement: - 
Local Independence: - 
Smooth: 4- 
Adaptive: + 

Minimize slivers: + 
Sharp features: + 

Thin features : + 


Image credit: http : // www.duang le,c^^^ 

http://www.oqre3d.org/tikiwiki/tiki- 

|ndex,php?pa9e=VolumeComponent 



L. 
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Cubical Marching Squares 

https://www.csie.ntu.edu.tw/~cyy/publications/papers/Ho2005CMS.pdf 

• Build an octree with error subdivision (similar to DMC) 

■ For any voxel (at any octree level): 

■ Unfold the voxel, look at each side individually 

■ Create a curve using marching cubes, tessellating if error > e 

■ Fold sides back together and triangulate 



Image credit: 

https://www.csie.ntu.edu.tw/~cyy/publlcations/ 

papers/Ho2005CMS.pdf 


Cubical marching squares treats each side as an independent 
problem, subdivides to match adjacent octree nodes or fit 
geometry better. 

Triangulation can use any method once you reassemble the 
box and extract the curves 
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Cubical Marching Squares 


Easy to Implement: - 
Local Independence: + 
Smooth: 4- 
Adaptive: + 

Minimize slivers: + 
Sharp features: + 

Thin features : + 


Image credit: htt ps : //upvoid.com / 

https://community.upvoid.com/uploacls/default/76/f220 

5 d4 j9 d ip9bbc8 . jpg 
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Windborne - Overview 

F(x,y,z) = composited, axis-aligned density grid 
• Each chunk has the list of features overlapping it 
■ On calculation, accumulate to chunks with boolean ops 



Windborne was a voxel game we worked on, mostly in 2013- 
2014. 

Chunks are just regions of space (16x16x16 for instance) that 
help manage vertex buffers, visibility, and streaming. 

Basic operation is "Here's a chunk, give me back the density 
grid (and materials)" 



Voxel game, regular grid because we found it was much easier 
for users to understand and interact with. It's interesting as 
well that you need a lot more detail and variation when you go 
smooth terrain; we initially started with something like MC 
source and it was very boring once the high frequency of 
individual voxels was removed. 
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Lessons - Compositing 


- Each feature is a layer 

■ Each layer either 
subtracts or overlays 

■ Alpha blend densities 

a n + a 6 (l - a a ) 

Image credit: 

https://en.wikipedia.org/wiki/Alpha_compositing 
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This was what worked best. Min/max/etc didn't for various 
reasons (dealing with surfaces that were almost aligned) 
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Lessons - Compositing 



■ Caves are a 
'subtract' layer 



Llllk 


29 



Lessons - Uniform representation 


■ Player wants objects + terrain to 
feel the same 

■ Collocate density grid and hand 
crafted voxel geometry 

■ Needs independence 

■ Draw sides of the voxel 




1 1 n 

Storing solid sides with our MC solution was pretty simple, and 
it meant we could draw them when there was a gap (to avoid 
seams), while ignoring them if the densities of neighboring 
voxels lined up 



Lessons - Uniform representation 


■ Represent density as 24 bit mask (3 bits for each 
corner), fits in with non-terrain voxels 

■ When well structured, VERY fast query operations 

Block: :Terrain IsSideFullySolid( BlockFace: :Face face ) const 
{ 

return (( Terrain.mVertexDensity & sBHsolidMaskForFace[ face ] ) *= sBHsolidMaskForFace[ face ] ); 

> 

] Block: :Terrain_IsBlockFullySolid() const 
{ 

return (( Terrain.mVertexDensity & 044444444 ) -* 0444444 4 4 ); 

} 


LI 


Note that for each voxel, we store the densities of all the 
corners. This is redundant, but important from a user 
perspective - independence was VERY important - don't edit 
one thing and then have something neighboring change. 
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Lessons - Uniform representation 



LIN 


Raise terrain under features to line up density field / non- 
density-field features 
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Lessons - Contour 


http://josiahmanson.com/research/contour indicator/ 


Image credit: 

http://josiahmanson.com/research/contour_lndicator/ 


Indicator 


llli 

MC underrepresents volumes - for instance if a single corner 
has a density of '1' and the other 7 have a density of 0, you'd 
expect about an eighth of the voxel to be filled in. Instead of 
that though, the traditional algorithm basically slices off a 
corner half way down three of the edges, whose volume is 
l/6 th of what it should be. 


With some simple math we can correct for this, and it 
mitigates the ridges that you will otherwise see. 
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Lessons - Density Grid 




■ Dynamic Lighting 
■ Light flow 
- AO 



Grids are very good for dynamic lighting- many modern 
lighting techniques simulate light flow over a grid, and we 
additionally have densities along each edge and materials. 
For our placed objects, we actually made a uniform model of 
density and material to unify lighting across everything. 
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Lessons - Density Grid 


■ Occlusion mapping 

(Heirarchial) 



Because we had a well structured grid, we built an octree for 
occlusion queries- we could keep track of whether everything 
in the octree node was solid or transparent, and make some 
simplifying assumptions for rendering as a result. 
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Lessons - Expose Parameters 
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Exposing the parameters of the voxelization algorithm allowed 
for some great effects to minimize some of the poor features 
of the algorithm (Marching cubes lack of high frequency 
features for instance). We let the normal smoothing 
parameter vary across materials and it really helps with 
varying the character of what you see. 




Lessons - Leverage Idiosyncrasies 



Remember the '2 voxel' nyquist frequency- well we had 
redundancy which gave us the ability to represent higher 
frequencies, if we used it. We were able to build terracing into 
our generation algorithm. 
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Questions? 


Email: michael@hiddenpath.com 


